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SYMBOLS 


B Tunnel width 

b Ratio of model chord to tunnel width (< C/B ) 

C Model chord 

Cd Sectional drag coefficient 

C p Coefficient of pressure 

c s Source/sink spacing distance 

D Sectional drag 

Q Source strength 

q Tunnel dynamic pressure (|pu 2 ) 

Uqq Freestream velocity 

A u Perturbation velocity 

Ax Width at half height of symmetric signature 
x Axial spatial coordinate 

y Lateral spatial coordinate 

e Centerline interference velocity; proportionality constant 

p Density 

Subscripts: 

B Horizontal bouyancy 

c Corrected for blockage effects 

m Measured 

p Peak 

s Body /bubble (symmetric) 

w Wake (antisymmetric) 

0 Origin 

oo Freestream 


in 



SUMMARY 


The Wall-Pressure Signature Method for correcting low speed wind tunnel data to free-air 
conditions has been revised and improved for two-dimensional tests of bluff bodies. The method 
uses experimentally measured tunnel wall pressures to approximately reconstruct the flow field 
about the body with potential sources and sinks. With the use of these sources and sinks, the 
measured drag and tunnel dynamic pressure are corrected for blockage effects. Good agreement is 
obtained with simpler methods for cases in which the blockage corrections were about 10% of the 
nominal drag values. 


INTRODUCTION 

In recent wind tunnel tests of the downloads on the wings of the XV-15 Tilt-Rotor during 
take off and hover, large blockage effects were encountered (ref. 1). The validity of conventional 
theoretical blockage corrections used in the data reduction of these tests has not been established. 
As a result it was decided to use an alternative correction method, based more on experimental 
information, for a second series of tests. The chosen method, called the wall-pressure signature 
method, uses pressure distributions on the tunnel walls measured during the wind tunnel tests to 
better predict blockage corrections. 

In the wall-pressure signature method the flow field about the body is approximated using the 
superposition of flows associated with a set of sources and sinks. The strengths and positions of 
these sources and sinks are determined so as to reconstruct the measured velocity distribution on 
the tunnel walls. Once determined the effect of the tunnel walls on the measured drag and dynamic 
pressure at the model is estimated and appropriate blockage corrections made. 

The method was originally devised for three-dimensional (3-D) tunnel setups by Hackett, Wils- 
den, and Lilley (ref. 2) and has been revised for the two-dimensional (2-D) wind tunnel tests of 
the XV-15 wings. 

The present report outlines the 2-D version of the wall-pressure signature method and gives a 
comparison between the blockage corrections obtained using this method and those of reference 1. 
In addition, descriptions of the programs used in the calculation of the blockage corrections are 
given in appendicies A, B and C. 
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WALL-PRESSURE SIGNATURE METHOD 


Figure 1 shows a typical experimental setup for a 2-D wind tunnel. During the wind tunnel 
tests, the pressure distribution along the tunnel walls is one of the measurements recorded. This 
pressure distribution is converted to incremental or perturbation velocities about the freestream 
velocity (Uoo) by use of the definition of dynamic pressure 

^ = y'l-AC p -l ( 1 ) 

where A C p is the net pressure coefficient after the wall pressure coefficients for the empty tunnel 
have been subtracted off. 

The resulting incremental velocity distribution is assumed to consist of the superposition of the 
velocities for two flow fields, one symmetric and the other antisymmetric (fig. 2). The symmetric 
signature, modeling the body and its separation bubble, is constructed from a point source/sink 
pair (±Q») at a distance c, apart. The antisymmetric signature, modeling the viscous wake of the 
body, is obtained from a single point source (Q w ) located at the peak of the symmetric velocity 
distribution. 


I. Antisymmetric Signature Modeling (Wake) 

For a point source of strength Q located at (xo, yo), the x-component of induced velocity at an 
arbitrary location (x, y) is given by 

_ Q r x-xo 1 r „v 

Au 4k [(x — x 0 ) 2 + (y — yo) 2 J * K ) 

Although the wake signature is modeled by a single source Q w , a sink of equal strength at some 
downstream location must accompany it to ensure mass conservation. Further, the effects of the 
walls on the flow field can be simulated by the superposition of an infinite row of image systems 
as shown in figure 3. Thus, the velocity increment on one of the walls (y = ±B/ 2), resulting from 
the wake source/sink pair, is 

{ i f q \ r ^ — X2 x "*• £5 

ST \U^b) [(x-x 2 )2 + („+ 1/2)2 - (x-x 5 )2 + (n+ 1/2)2 

where the source is located at (x 2 , 0) and the sink is located at (xg, 0). The barred distances are 
nondimensionalized by the tunnel width B. 

This equation is too cumbersome to calculate at each wall port location, especially considering 
that equation 3 is very slowly convergent. Hence, it will be approximated by a hyperbolic tangent 
function 

(jj—'j = Aijl + tanh[A 2 (x-x 2 )] j (4) 
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where the constants A\ and Ai are determined from numerical analysis of equation 3. Note that 
the downstream asymptote is 2A\ and the slope at x = £2 is A 1 A 2 . 

The downstream asymptote, taken to be the peak velocity at x — 1/2 (£2 + £ 5 )) and the slope 
at X 2 of equation 3 vary only slightly for the downstream sink located in the range 10 + 22 < £5 < 
1000 + X 2 - Using a large number of image systems (» 50,000), the summation terms (in braces) 
approach n and the slope at x = X 2 is 4.800. Thus the constants A\ and A 2 become 

Ai =Kife)' A2=3 - 056 - (5) 

In addition, since A\ is half the asymptotic downstream velocity, the wake source strength is 
then given by 



Note that equation 6 is in conflict with the results of Hackett (ref. 2) by the factor of 2. The 
effect of this discrepency between the present analysis and that of Hackett’s on the final blockage 
corrections is further detailed in section II of Results. Because of this discrepency a more detailed 
description of this analysis is also given in appendix D. 


II. Symmetric Signature Modeling (Body /Bubble) 


Once the wake signature is determined, it is subtracted from the measured wall velocities, 
leaving that portion due to the body /bubble. The resulting symmetric signature is curve fitted by 
a parabola 

p) 

From this curve fit, the peak velocity and position are determined along with the width at half 
height. The data is filtered such that the points used for the curve fit properly model the upper 
half of the peak. From equation 7 these are given by 


Peak Velocity : ( Au/Uco) max = 7 - ft 2 / 4a 
Peak Position : x p = — /3/2a 


Width at Half Height : 


- =2 > ( AU/M, 


2a 


( 8 ) 

( 9 ) 

( 10 ) 


Once the parabola is determined, it becomes an inverse problem to find the source/sink strength 
and spacing which corresponds to that distribution. This task is accomplished by using x as input 
for two interpolation tables (appendix C). The output of these tables is the source/sink spacing 
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(c,) and the maximum velocity normalized by source strength (AuB/Q,) max , which is used with 
equation 8 to obtain the symmetric source/sink strength 


Q» _ (Au/U m ) max 
UooB (Au B/Q,) max ’ 


The source and sink positions are given by 

1 1 
x$ = Xp- -c,, £4 = £p+ -c,. 


(11) 


( 12 ) 


« 

HI. Iteration Procedure 

Since the position of the wake source (£ 2 ) is not known prior to the use of equation 4, a value 
must first be assumed and then iterated upon. 

Initially, the wake source position is assumed to be at the model location. The downstream 
asymptote is then determined from the data (see section I of Results), and the antisymmetric 
velocity distribution is calculated from equation 4. The result is subtracted off of the measured 
velocity distribution, and the resulting symmetric signature is curve fitted by the inverse parabola. 
If the peak position (eq. 9) is not sufficiently close to the assumed wake source position, a new 
value of x p is chosen and the procedure is repeated until convergence is obtained. 

Once this process is complete, the symmetric source/sink strengths and positions are obtained 
as outlined in section II, and the wake source strength is calculated by equation 6. 

If the curve fit yields a divergent result, the symmetric signature is smoothed by replacing the 
value at each point by the average of the point and its two immediate neighbors. This averaging is 
done only once. If convergence is not obtained after the smoothing operation, the program defaults 
to a lower order linear theory, known as HensePs method (ref. 3), to calculate the centerline 
interference velocity due to the symmetric signature. 


IV. Centerline Interference Velocities 

Once the source strengths and positions are known, the total centerline interference velocity 
distribution is obtained by superposing the effects of each of the three sources. For each source 
the interference velocity is found by using the position and strength as inputs for an interpolation 
table (appendix C). 

As stated previously, if convergence cannot be obtained, the antisymmetric interference velocity 
is calculated as normal, but the symmetric signature contribution is obtained by HensePs method 
(ref. 3). For 2-D source flow in a channel, HensePs result is that for x-positions relatively close to 
the source position, the ratio of wall velocity to centerline interference velocity is simply 3. 
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V. Blockage Corrections for Dynamic Pressure and Drag 


Given the interference velocity distribution along the tunnel centerline, the maximum velocity 
((max) is found and used to correct the tunnel dynamic pressure 

9c = 9m(l + Cm«) 2 (13) 

where q m is the measured tunnel dynamic pressure and q e is the corrected dynamic pressure. 
Corrections to the measured drag coefficient include both effects of the tunnel q corrections and 
horizontal bouyancy. 

In the original program of Hackett, horizontal bouyancy drag corrections were obtained by a 
method which uses the axial C p gradient at the body. However, in the present test cases the axial 
extent of the model is not precisely known because of the high angles of attack. As a result it was 
decided that a second method should be used. This method, called the “ puQ ” method, recognizes 
from a momentum balance that the bouyancy drag on the model (A Db) is given by 

A D b = p{- Au+ Q, + AuJQ t - Au+<5 W ) (14) 

where Au* is the induced velocity of each source. However, as mentioned before, a downstream 
wake sink is necessary for continuity; the drag of all four sources is then zero: 

ADtot = p(-Au+Q, + AuJQ, - A u+Q w + A u~Q w ) = 0. (15) 

Thus, taking the difference between equations 14 and 15 gives 

A D b = -pAu~Q w . (16) 


Now, A u w is half the asymptotic velocity of the wake source which from equation 6 becomes 



Then the drag coefficient correction for horizontal bouyancy is given by 



where C is the model chord. 

The total drag correction including dynamic pressure and bouyancy effects is 

Cd c = {C Dm + AC DB )(q m /q e ) (19) 

where Cjj m is the measured drag coefficient and Cd c is the final, corrected drag coefficient. 
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VI. Reconstruction of Wall Velocities 


Among the output from the program (appendix A) is a graphical comparison of the input tunnel 
wall velocities to those calculated using the source/sink strengths and positions. This gives the 
user an indication of the accuracy that the tunnel wall velocities are reconstructed by the pressure- 
signature method. These reconstructed wall velocities are calculated in the following manner. Like 
the wake source (eq. 4) the effect of the symmetric source/sink pair is modeled by two hyperbolic 
tangent functions. The resulting reconstructed wall velocity signature is then given by 


^ = K^) i 1 + ■ * 2)| } 

+ ^ { tanh [^ 2 (* - * 3 )] - tanh[A 2 (* - * 4 )] }. 


( 20 ) 


I 

DISCUSSION OF RESULTS 


I. Downstream Asymptotic Velocity 

To calculate the antisymmetric velocity distribution, the asymptotic velocity is required. How- 
ever, this value is not always accurately known because of such things as tunnel length restrictions 
and data spread. Thus, the sensitivity of the results to the choice of this asymptotic velocity must 
be investigated. 

Figure 4 shows the measured pressure distribution for a test of the 30-cm triangle with apex 
forward (ref. l). This case was analysed several times using different ports for the asymptotic 
velocity. The results, shown in table I, are typical of tests in this series. 

As is shown in the table, the effect of the asymptote on the final dynamic pressure and drag 
corrections is small, even though there is a relatively large spread in the parameters of the individual 
sources (especially in their positions). 

Using the search option of the program, the results would correspond to port 14 since it has 
the lowest velocity of the last few ports. Note also that the current version of the program will not 
use the last port as asymptote if the second to last port has a lower velocity. 


II. Comparison with Hackett’s Corrections 

To see the effect of the discrepancy in equation 6 between Hackett’s method and the present 
analysis on the final corrections, the aforementioned case was rerun using Hackett’s version of 
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Table I: Comparison of Results Using Different Asymptotes 


Port 

%Aq a 


wsm 

W33SM 

Q»/ub 

Q W /UB 

x 2 

C 8 

mm 

5.60 

hee 1 

0.0275 

-0.0923 

0.0265 

Mjgm 




5.69 

MM 


-0.0969 

0.0259 

EH 


0.531 


5.79 

-11.6 


-0.1014 

0.0232 

0.142 


0.535 

14 

5.41 

-10.2 

0.0265 

-0.0855 

0.0291 

0.131 

0.065 

0.715 

percent change in 

dynamic pressure 

- 7m)/ 



^percent change in drag coefficient (Cp c — Cj) m )fCv m 
‘solution on verge of divergence 


equation 6. Table II shows the comparison between the corrections for the two methods for the 
last entry in table I. 


Table II: Comparison of Corrections for the Two Versions of Equation 6 


Equation (6) 
Coefficient 

%Aq 

%AC d 

f max 

AC Db 

Q'/ub 

QJUB 

X 2 

m 

2.0 (present) 

■ 

-10.21 




0.1307 



1.0 (Hackett) 


-5.28 

Iwmll 


0 0291 

0.0653 

0.065 

\ 0.716 I 


Table II shows that Hackett’s corrections are approximately half of those using the present 
analysis of the wake strength to asymptotic velocity relationship. 


Ill* Comparison with Conventional Blockage Corrections 

As mentioned earlier, the reason for the present investigation was the uncertain validity of the 
conventional corrections used in reference 1 for the unusually large blockage effects found there. 
Here those corrections are compared with the present method for the triangle case of section I. In 
reference 1, the blockage correction formulas used were 

?C = ?m( l + e&), C D C = C Dm (l - ebC Dm ) (21) 

where b is the ratio of model width to tunnel width, and e was estimated to be 0.65 ± 0.05. For 
the triangle case of section I, b = 0.10 and the measured drag coefficient was Cr> m — 1.582. The 
corrections using these values are compared with those of the present method in table III. 

The close comparison of the two methods shown in table III gives an indication that the cor- 
rections used in reference 1 were, in fact, valid despite the magnitude of the blockage effects. 
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Table III: Comparison with Convetional Blockage Corrections 


Method 

%Ag 

% ac d 

Pressure-Signature 
McCroskey (ref. 1) 

difference 

5.41 

6.50 

1.11 

-10.21 

-10.28 

0.07 


CONCLUSIONS 

The wall-pressure signature method for correcting low speed wind tunnel data to free-air condi- 
tions has been revised and improved for 2-D tests of bluff bodies. The method uses superposition 
of the flow fields associated with a set of three linear potential sources to approximate the flow 
about the body in the presence of the wind tunnel walls. Strengths and positions of the sources 
is determined so as to reconstruct the velocity distribution on the tunnel walls, which is obtained 
from measured pressure distributions taken during the wind tunnel tests of the model. With the 
use of these sources and sinks, the measured drag and tunnel dynamic pressure are then corrected 
for blockage effects. 

This method has been used to apply blockage corrections of 2-D wind tunnel tests performed 
on the downloads on the wings of the XV-15 Tilt-Rotor during take off and hover. In these tests 
the blockage corrections were on the order of 10% of the measured drag values. The corrections 
obtained with this method were found to be in good agreement with the simpler methods used in 
reference 1. 
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APPENDIX A 


MAIN PROGRAM DESCRIPTION: BLKAGE2D 

The main blockage calculation program (BLKAGE2D) is a revised version of the 3-D code 
supplied by Hackett, Wilsden, and Lilley (ref. 2) for the CDC 7600 at NASA Ames Research 
Center. 

The main operational difference between the two codes, other than conversion from 3-D to 2-D, 
is that if the symmetric signature cannot successfully be curve fitted, its points are smoothed. Then 
another curve fit is tried and if again unsuccessful, the program branches to Subroutine PUNT. In 
PUNT the wake, or antisymmetric portion, is kept and the symmetric portion is modeled by the 
Hensel computation. The impetus behind these changes is that it was found that for the present 
experiment, the wake signature alone fit the data well, leaving little more than experimental scatter 
for the symmetric portion in many cases. 


Inputs 


The program needs two separate inputs; the first is the individual run inputs and the second is 


the lookup charts. 


Run Input (Sxxx.TMP) 


Individual run input can again be broken down into two parts: main input, which is inputted 
once per run, and frame input for each frame of data. These two inputs have the following form. 
Unless otherwise specified, all input is in free format: 


Main Input: 


I I) RUNUM (A3) 

2) BTUN, CMOD, XMOD 

3) NWST, NBST 

4) XWST(I) 1=1, NWST 

5) XBST(I) 1=1, NBST 

6) LU, IUSES, IUSEE, ILIST, IDEBUG, ITAB 

7) IOPT 
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RUNUM = run number (3 digits) 

BTUN = width of tunnel fft) 

CMOD = chord of model (ft) 

XMOD = axial position of model (ft) 

NWST = number of wall pressure ports 
NBST = number of body pressure ports 
XWST(I) = axial position of Ith wall port 
XBST(I) = axial position of Ith body port 
LU = plotter device number (=0) 

IUSES = forward velocity asymptote 
= 0 - zero asymptote 

> 0 - velocity at IUSES port 
IUSEE = aft velocity asymptote 

= 0 - searches for smallest velocity after peak 
< 0 - velocity at IUSEE port from end 

> 0 - averages last IUSEE ports 
ILIST = additional output option 

= 0 - no added output 

<> 0 - distripution of C p and velocity along walls and C-L 
IDEBUG = debugging output (no output if = 0) 

ITAB = lookup charts output (no output if = 0) 

IOPT = next input option 

= 1 - new main input 
= 0 - new frame input 
=— 1 - end 


Frame Input: 

I I) ALPHA 

2) CPWST(I) 1=1, NWST 

3) CPEM(I) 1=1, NWST 

4) CPB(I) 1=1, NBST 

5) QU, PU 

6) CHDAT, CHTIM, CHILE, CHITE, CHVAR (2A8,3A15) 

7) IFRAME, CMUU, CLU, CDU, CMU 

8) IOPT 

ALPHA = angle of attack (deg) 

CPWST(I) = measured wall C p at Ith port 
CPEM(I) = empty tunnel wall C p at Ith port 
CPB(I) = measured body C p at Ith body port 
QU = measured dynamic pressure (psf) 

PU = measured static pressure (psfa) 

CHDAT = date of experiment (xx/xx/xx) 

CHTIM = time of experiment (xx:xx:xx) 

CHILE = description of leading edge 
CHITE = description of trailing edge 
CHVAR = variation description 
IFRAME = frame number 

CMUU, CDU, CLU, CMU = measured force coefficients (power, drag, lift, moment) 
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Charts Input (LOOKUP.TAB) 
Lookup charts input is in the same form as output from LOOKUP: 

I I) NDX, NX 

2) XDXOB(I) 1=1, NDX 

3) CSOB(I) 1=1, NDX 

4) XUFM(I) 1=1, NDX 

5) XXOB(I) 1=1, NX 

6) XUF(I) 1=1, NX 

7) AT(I) 1=1,3; AH(I) 1=1,2 

NDX = number of points in Charts I and II 
NX = number of points in Chart III 
XDXOB = width at half height 
CSOB = source/sink spacing 

XUFM = max velocity normalized by source strength 

XXOB = axial position in tunnel 

XUF = centerline interference velocity 

AT = tanh constants (A,- in equations 4 and 20) 

AT(1) = 3.056, AT(2)=AT(3)=0 
AH = Hensel constants 

AH(1) = 1/3, AH(2)=0 


Output (Sxxx.OUT) 

The program prints a main output (including a plot of measured and calculated wall velocities), 
along with three optional outputs depending on the values of ILIST, IDEBUG, and ITAB. In 
addition, a summary output is printed (Sxxx.SUM). 

The following is a description of the output variables for each. 
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Variable 


Output 


Description 


Main Output: 


EPS (MOD) 

EPSMOD 

EPS (MAX) 

EPSMAX 

X(MOD)/B 

XMOD 

X(MAX)/B 

XP 

XV/B 

XVOB 

BS/B 

BSOB 

DX/B 

DXOB 

CS/B 

CSOB1 

QS/UB 

QFS 

QW/UB 

QFW 

DCDWB 

DCDW 

US(MAX)/U 

UOUMAX 

UFM 

UFMAX 

H FACTOR 

HENSEL 

A5 

A5 

A6 

A6, AT(1) 

A7 

A7 

Additional Output: 

X/B 

XWST 

CP 

CPWST-CPEM 

u/u 

UOU 

UA/U 

UAOU 

US/U 

USOU 

uw/u 

uwou 

UP/U 

UPOU 

uv/u 

uvou 

EPS 

SIGUOU 

Lookup Chart Output : 


C-L interference velocity at model 
maximum C-L interference velocity (e m ax) 
axial position of model 

axial position of peak symmetric' velocity (£ 2 ) 
axial position of wake source (Q w ) 
not used in 2-D 
width at half height 
source/sink spacing 
solid body source strength (Q t ) 
wake source strength (Q w ) 
bouyancy drag correction 
maximum symmetric velocity on wall 
maximum source-strength-normalized velocity 
ratio of peak symmetric wall velocity to 
maximum C-L interference velocity 
half of asymptotic downstream velocity 
tanh constant (A, in eqs. 4 and 20) 
not used in 2-D 


axial position of wall ports 
measured zeroed C p along wall 
measured incremented velocity (plotted) 
antisymmetric (wake) velocity along wall 
symmetric (body /bubble) velocity along wall! 
computed wall velocity (plotted)' 
computed C-L interference velocity 
computed C-L interference velocity 
computed C-L interference velocity 


(body! 
(wake) 
(total) 


- same as input format 


Debug Output: 


— see code (Subroutine EPSCAL) 



APPENDIX B 


AUXILIARY PROGRAMS DESCRIPTION 

Two additional programs are needed for running BLKAGE2D. The first is a preprocessor, 
BLSETUP, and the second is a routine to setup the interpolation table for the empty tunnel wall 
pressure distribution. 


Program BLSETUP 

The main purpose of BLSETUP is to reduce wind tunnel data into pressure and force coefficients 
and then output the reduced data in a form which can be read in by BLKAGE2D. 

The program is set up so that a subroutine reads in and reduces individual frames. This 
subroutine is taken from an off-line analysis routine called NEWA: 

Inpifts: Sxxx.DAT - wind tunnel data 

Syyy.EMP - empty tunnel C p interpolation table 

Outputs: Sxxx.TMP - BLKAGE2D input 


Program BLEMPT 

Because the empty tunnel pressure distribution may change with tunnel dynamic pressure (i.e., 
Reynolds effects), Program BLKAGE2D was changed so that empty tunnel C p s are read in for 
individual frames rather than once per run. Thus, BLSETUP must output empty tunnel C p s 
which are appropriate for each frame’s dynamic pressure. This task is accomplished by the use of 
an interpolation table in which tunnel q is the independent variable. 

Program BLEMPT constructs this table. Since empty tunnel data is stored in the same format 
as normal runs, Program NEWA is again uses for data input and reduction. 

After all frames of the empty tunnel run are input and reduced, the program prompts for which 
frames to average. 

Inputs: Syyy.DAT - empty tunnel data 

Outputs: Syyy.EMP - empty tunnel C p interpolation table 

It should be noted that the current versions of BLSETUP and BLEMPT disregard pressures 
from port 12 and replace them with the average of ports 11 and 13. 
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APPENDIX C 


LOOKUP CHARTS 

After the antisymmetric signature is subtracted from the measured wall velocity distribution, 
the resulting symmetric signature is curve fit by an inverted parabola. Then an inversion process 
is performed to obtain the source/sink strengths and positions which correspond to this parabolic 
distribution. This process is accomplished using interpolation tables (Charts I and II). 

Further, after all source strengths and positions have been found, the centerline interference 
velocity must be found. This process also is accomplished by the use of an interpolation table 
(Chart III). 


Charts I and II 


Using the width at half height (Ax) as input, Chart I outputs the source/sink spacing (c,), and 
Chart II outputs the maximum velocity normalized by the source strength (A uB/Q,) max . 

Considering a source/sink pair located at (0,0) and (c„0) in a tunnel (fig. 2), where the tunnel 
walls are simulated by a singly infinite row of image systems, the incremental velocity at a location 
on one of the walls (y = ±5/2) is given by 

Au = J_ ( Q> \ y? [ x x-c, 1 

Uoo 4x \UooBj n ~' 00 [ x 2 + (n + 1/2) 2 (x - c,) 2 + (n + 1/2) 2 J 


or 


A «5 

Q» 


_i_ r x x-c, 

2x “ [x 2 + (n + 1/2) 2 (x - c,) 2 + (n + 1/2) 2 


where barred distances are normalized by tunnel width B . 


(Cl) 


Charts I and II are constructed by the following procedure: For a given range of c sy the 

maximum source-strength-normalized incremental velocity (A uB/Q 9 ) is determined by evaluating 
equation Cl at the midpoint x = \c 9 . The width at half height is determined by iteration. This 
is done by evaluating equation Cl for different values of x until the position x = X 1/2 at which 
A uB/Q 9 = l/2(AuB/<2,) max is found. Then the width at hair height is given by 


Ax = 2(x 1/2 - 



(C2) 
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Chart IH 


Using x position as input, Chart III gives the centerline interference velocity (normalized by 
source strength) caused by the presence of tunnel walls on a single source of strength Q located on 
the centerline. 

As in equation Cl, the wall effects are simulated by a singly infinite row of image systems. 
However, in this case the equation is evaluated on the centerline (y = 0) and only the effects of the 
image systems are included. 



C-L 


J_ ( Q \ £5 [ x-xo 1 
2jt UooJB/^L(*-*o)* + n a J* 


Chart III is constructed by evaluating equation C3 over a given range of x positions. 


(CS) 


Program Description (LOOKUP) 

The program used for construction of the lookup charts is called LOOKUP. Its inputs are as 
follows: 


Charts I and II 
Chart III 


(1) minimum and maximum values of c 8 

(2) number of points in Chart I and II 

(3) minimum and maximum values of x 

(4) number of points in Chut III 

(5) number of image systems used in calculation 

(6) iteration error parameter 


It is suggested that a minimum value for c 8 of not less than 0.05 be used. Also, since equations 
Cl and C3 are slowly convergent series, the number of image systems should be on the order of 
10 5 or 10 6 . 
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APPENDIX D 


ANALYSIS OF WAKE VELOCITY APPROXIMATION 


As stated in the text, there is a disagreement between Hackett’s results and the present analysis 
in the relationship between the asymptotic velocity and the strength of the wake source. Therefore, 
it has been decided to detail the analysis of the wake velocity distribution and its approximation 
by a hyperbolic tangent function. 

As stated in equation 3, the induced velocity distribution on the tunnel wall caused by the wake 
source/sink pair is given by 

2 f _ 1 _ / Qw \ y? r X-X2 x-x 5 ' 

(47T \UooBj £? q [(x — X2) 2 + (» + 1/2) 2 (x - x 5 ) 2 + (n+ 1/2) 2 . 

where the factor of 2 comes from the fact that the summation terms account for the image systems 
on one side of the wall only (the influence of the images is symmetric about the wall). 

Numerical experiments were run on the summation terms alone assuming the wake source to 
be at X2 = 0. Table D-I shows the results for the downstream asymptote and slope at x = X2 for 
various sink locations (£5) and number of image systems used. The downstream asymptote is taken 
to be the peak velocity at the midpoint between the source and sink (x == |xs). 




Table D-I: Summation Terms for Asymptote and Slope 


X 5 

Number of 
images 

Asymptote 

Slope @ X2 

10 

1,000 

3.132 

4.801 

10 

10,000 

3.140 

4.801 

10 

100,000 

3.140 

4.800 

100 

10,000 

3.132 

4.801 

100 

100,000 

3.137 

4.798 

1000 

10,000 

3.133 

4.801 

1000 

100,000 

3.128 

4.804 


Table D-I shows that the asymptote and slope are nearly constant over a large range of sink 
locations and number of images used. From the table the asymptote can be taken as 7r and the 
slope at X2 as 4.800. This gives the asymptotic velocity 



(D2) 
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(D3) 


and the slope at the source location 

±(*±)\ = J-(^[ 4 . 800 l 

dx \UooJ |, 3 2t \UooBj [ J 


2A ( Qvj \ 

* \UooBJ' 


The actual velocity distribution (eq. Dl) is approximated by a hyperbolic tangent function as 


in equation 4. 

= Ai|l + tanh[A 2 (s-* 2 )]}- 

Noting that tanh(x) — ► 1 as x — ► oo, the asymptote of equation D4 is 

(D4) 


HI 

l 

8 

II 

to 

U- 

(D5) 

The slope at x = x% is given by 


d / Au\ 1 

di\uZ)\ 

= A 1 A 2 s 1 + tanh [A 2 {x — 12 )] \ = A 1 A 2 . 

1 2 l J Jt=Z 3 

(D6) 

Thus, from equations D2 and D5, A\ is given by 



*=K«S>) 

(D7) 

and from equations D3, D6, and D7, A 2 becomes 



A 2 = 3.056 

(D8) 
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PRESSURE PORTS 



Figure 1. - Sketch of 2-D wing model in tunnel with wall pressure ports 
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Figure 2. - Approximation to the flow field in the wind tunnel using potential sources and sinks 
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Figure 3. - Representation of wind tunnel walls by image systems 
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Figure 4. — Pressure distribution along wind tunnel walls for a model with triangular cross 

section; apex forward, b = 0.10 
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